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We report upon the numerical computation of the Eu- 
ler characteristic x (a- topologic invariant) of the equipoten- 
tial hypersurfaces E„ of the configuration space of the two- 
dimensional lattice if"^ model. The pattern xC^v) vs v (po- 
tential energy) reveals that a major topology change in the 
family {E„}„gR is at the origin of the phase transition in the 
model considered. The direct evidence given here - of the rel- 
evance of topology for phase transitions - is obtained through 
a general method that can be applied to any other model. 
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Suitable topology changes of equipotential submani- 
folds of configuration space can entail thermodynamic 
phase transitions. This is the novel result of the present 
Letter. The method we use, though applied here to a par- 
ticular model, is of general validity and it is of prospective 
interest to the study of phase transitions in those systems 
that challenge the conventional approaches, as it might 
be the case of finite systems (like atomic and molecular 
clusters), of off-lattice polymers and proteins, of glasses 
and in general of amorphous and disordered materials. 
Let us begin by giving a theoretical argument and then 
proceed by numerically proving its truth for the 2d lat- 
tice ip'^ model. Consider classical many particle systems 
described by standard Hamiltonians 
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where the (p, q) = (j>i, . . . ,pn, qi, ■ ■ ■ , qN) coordinates 
assume continuous values and V{q) is bounded be- 
low. The statistical behaviour of physical systems de- 
scribed by Hamiltonians as in Eq.(|l|) is encompassed, in 
the canonical ensemble, by the partition function in phase 
space 
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term is written using a co-area for- 
V labels the equipotential hypersur- 



faces St, of configuration space, — {(qi, . . . , gjv) G 
R^|T^(gi, . . . , gAf) = v}. Equation (^) shows that for 
Hamiltonians (|l|) the relevant statistical information is 
contained in the canonical configurational partition func- 
tion Z'^ = J Hdgi exp[— /3F((7)]. Remarkably, Z'^ is de- 
composed - in the last term of Eq.(^ - into an infinite 
summation of geometric integrals, da /|| Vy|| , defined 
on the {Si,}^gR. Once the microscopic interaction poten- 
tial V{q) is given, the configuration space of the system 
is automatically foliated into the family {Et,}t,gR of these 
equipotential hypersurfaces. Now, from standard statis- 
tical mechanical arguments we know that, at any given 
value of the inverse temperature /?, the larger the num- 
ber N of particles the closer to = are the mi- 
crostates that significantly contribute to the averages - 
computed through Zjv(/3) - of thermodynamic observ- 
ables. The hypersurface S^i^ is the one associated with 
up = [Z'^)-'^ j Y{dqiV{q)e-l^^^'^\ the average potential 
energy computed at a given (3. Thus, at any /3, if is 
very large the effective support of the canonical measure 
shrinks very close to a single Ei, = . Hence, and 
on the basis of what we found in let us make ex- 

plicit the Topological Hypothesis: the basic origin of a 
phase transition lies in a suitable topology change of the 
{Su}, occurring at some Vc- This topology change induces 
the singular behavior of the thermodynamic observables 
at a phase transition. By change of topology we mean 
that {'£iv}v<v^ are not diffeomorphic to the {S]^}t,>t,^ 
In other words, the claim is that the canonical measure 
should "feel" a big and sudden change - if any - of the 
topology of the equipotential hypersurfaces of its under- 
lying support, the consequence being the appearence of 
the typical signals of a phase transition, i.e. almost sin- 
gular (at finite N) energy or temperature dependences of 
the averages of appropriate observables. The larger TV, 
the narrower is the effective support of the measure and 
hence the sharper can be the mentioned signals, until true 
singularities appear in the N oo limit. This point of 
view has the interesting consequence that - also at finite 
- in principle different mathematical objects, i.e. man- 
ifolds of different cohomology type, could be associated 
to different thermodynamical phases, whereas from the 
point of view of measure theory the only mathemati- 
cal property available to signal the appearence of a phase 
transition is the loss of analyticity of the grand-canonical 
and canonical averages, a fact which is compatible with 
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analytic statistical measures only in the mathematical 
— > oo limit. In order to prove or disprove the con- 
jectured role of topology, we have to explicitly work out 
adequate information about the topology of the members 
of the family {S„}^gR for some given physical system. 
Below it is shown how this goal is practically achieved by 
means of numerical computations. As it is conjectured 
that the counterpart of a phase transition is a breaking of 
diffeomorphicity among the surfaces , it is appropriate 
to choose a dijfeomorphism invariant to probe if and how 
the topology of the changes as a function of v. This is 
a very challenging task because we have to deal with high 
dimensional manifolds. Fortunately a topological invari- 
ant exists whose computation is feasible, yet demands 
a big effort. This is the Euler characteristic, a diffeo- 
morphism invariant, expressing fundamental topological 
information In order to make the reader acquainted 
with it, we remind that a way to analyze a geometrical 
object is to fragment it into other more familiar objects 
and then to examine how these pieces fit together. Take 
for example a surface E in the euclidean three dimen- 
sional space. Slice S into pieces that are curved triangles 
(this is called a triangulation of the surface) . Then count 
the number F of faces of the triangles, the number E 
of edges, and the number V of vertices on the tesselated 
surface. Now, no matter how we triangulate a compact 
surface S, x(^) = F — E + V will always equal a con- 
stant which is characteristic of the surface and which is 
invariant under diffeomorphisms : S — > S'. This is the 
Euler characteristic of E. At higher dimensions this can 
be again defined by using higher dimensional generaliza- 
tions of triangles (simplexes) and by defining the Euler 
characteristic of the n-dimensional manifold S to be 

n 

X{^) = ^(-l)''(#of "faces of dimension k"). (3) 

fe=0 

In differential topology a more standard definition of 
X(S) is 
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where also the numbers bk ~ the Betti numbers of S - are 
diffeomorphism invariants |^ . While it would be hopeless 
to try to practically compute x(S) from Eq.(P in the case 
of non-trivial physical models at large dimension, there 
is a possibility given by a powerful theorem, the Gauss- 
Bonnet-Hopf theorem, that relates x(S) with the total 
Gauss-Kronecker curvature of the manifold, i.e. |nO| 



x(^) =7 J Kg da 



(5) 



which is valid for even dimensional hypersurfaces of eu- 
clidean spaces [here dim(I]) — n = — 1], and 



where: 7 = 2/Vol{E>i) is twice the inverse of the vol- 
ume of an n-dimensional sphere of unit radius; Kg is 
the Gauss-Kronecker curvature of the manifold; da = 
y^det(g)dx^dx^ ■ ■ ■ dx" is the invariant volume measure of 
S and g is the Riemannian metric induced from M.^ . Let 
us briefly sketch the meaning and definition of the Gauss- 
Kronecker curvature. The study of the way in which an 
ri-surface E curves around in is measured by the way 
the normal direction changes as we move from point to 
point on the surface. The rate of change of the normal 
direction ^ at a point a; g E in direction v is described by 
the shape operator Lxiv) = — Vv^, where v is a tangent 
vector at x and Vv is the directional derivative, hence 
Lxi'v) = — (V^i • V, . . . , V^„+i • v); gradients and vectors 
are represented in M^. As is an operator of the tan- 
gent space at x into itself, there are n independent eigen- 
values Ki{x), . . . , Kn{x) which are called the princi- 
pal curvatures of S at x. Their product is the Gauss- 
Kronecker curvature: Kg{x) — Y['i=i'^ii^) ~ 'iet{Lx). 
The practical computation of Kg for the equipotential 
hypersurfaces proceeds as follows. Let |= Vy^/||Vy|| 
be the unit normal vector to at a given point x, and 
let {vi, . . . , v„} be any basis for the tangent space of E„ 
at X. Then [llll 
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Let us now consider the family of {E^}„gH associated 
with a particular physical system and show how things 
work in practice. We consider the so-called (p"^ model on 
a d-dimensional lattice Z'' with d = 1,2, described by the 
potential function 



V = 
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^ 2 



Jiq^-qk)^ (7) 



where (ik) stands for nearest-neighbor sites. This sys- 
tem has a discrete Z2-symmetry and short-range inter- 
actions; therefore, according to the Mermin- Wagner the- 
orem, in c? = 1 there is no phase transition whereas in 
d = 2 there is a symmetry-breaking transition of the same 
universality class of the 2d Ising model. Independently 
of any statistical measure, let us now probe, by com- 
puting x(^-u) V according to Eq.(^, if and how the 
topology of the hypersurfaces E„ varies with v. To this 
aim we first devised an algorithm of MonteCarlo type by 
constructing a Markov chain on any desired surface E^,. 
This is obtained by means of a "demon" algorithm cor- 
rected with a projection technique JT^ which provides a 
simple and efficient method to constrain a random walk 
on a level-hypersurface, here, of the potential function. 
Each new step so obtained on E^, represents a trial step 
which is accepted or rejected according to a Metropolis- 
like "importance sampling" criterion |1J] adapted to the 
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weight ^/deUjj). With any MonteCarlo scheme we can 
actuaUy compute densities, that is we can only estimate 
/j, Kcda/ /j, da, the average of Kq, rather than its 
total value ^ on S^,. Hence the need for an estimate 
of AreaCSy) — da as a function of v. To this aim 
we worked out a geometric formula that links the rela- 
tive variation of Area(Yiv) with respect to an arbitrary 
initial value AreaCEyg), to another MonteCarlo average 
on 5]„: {Mi/\\VV\\)l^c where Mi = ^ ELi is the 
mean curvature of E„ [|5j. Thus the final outcomes of 
our computations are the relative variations of the Eu- 
ler characteristic. The computation of Kq at any point 
X € E„ proceeds by working out an orthogonal basis for 
the tangent space at x, orthogonal to ^ = VV^/||Vy||, by 
means of a Gram-Schmidt orthogonalization procedure. 
Then Eq.(||) is used to compute Kg at x. On each 
we sampled 1 • 10^ — 3.5 • 10^ points where we computed 
Kg- This number of points was varied, and several initial 
conditions were also considered in order to check the sta- 
bility of the results. The computations were performed 
for dim(E„) = 48, 80 (i.e. iV = 7 x 7, 9 x 9) and with 
the choice A — 0.6, fi^ — 2, J ^ 1 for the parameters 
of the potential. In order to test the correctness of our 
numerical "protocol" to compute x{^v)j ^-^d to assess 
its degree of reliability, we checked the method against a 
simplified form of the potential V in Eq. (^, i.e. with 
A = J = 0, /i" = — 1. In this case the E^, are hyperspheres 
and therefore x(§") = 2 for any even n. Area(§") is an- 
alytically known as a function of the radius y/v, therefore 
the starting value Area{T,yg ) is known and in this case we 
can compute the actual values of xC^v) instead of their 
relative variations only. In Fig.l we report xi^v = §") vs 
v/N for iV = 5 X 5, the results are in agreement with the 
theoretical value within an error of few percents, a very 
good precision in view of the large variations of xC^v) 
that are found with the full expression (0) of V. In Fig.2 
we report the results for the Id lattice, which is known 
not to undergo any phase transition. Apart from some 
numerical noise - here enhanced by the more complicated 
topology of the E^ when A, J - a monotonously (in 
the average) decreasing pattern of x{u/N) is found. Since 
the variation of xi'^/^) signals a topology change of the 
{El,}, Fig. 2 tells that a "smoothly" varying topology 
is not sufficient for the appearence of a phase transition. 
In fact, when the 2d lattice is considered, the pattern 
of xi^/N) is very different: it displays a rather abrupt 
change of the topology variation rate with v/N at some 
Vc/N. This result is reported in Fig. 3 for a lattice of 
N ~ 7 X 7 sites, and in Fig. 4 for a larger lattice of 
= 9 X 9 sites The question is now whether the 

value Vc/N, at which x(w/iV) displays a cusp, has any- 
thing to do with the thermodynamic phase transition, i.e. 
we wonder if the effective support of the canonical mea- 
sure shrinks close to E„=:t,^ just at (3 = 1/Tc, the (inverse) 
critical temperature of the phase transition. The answer 



is in the affirmative. In fact, the numerical analysis in 
Refs. shows that - with A = 0.6, = 2, J = 1 - 

the function jf{V){T) and its derivative signal the phase 
transition at jf{V) ~ 3.75, a value in very good agree- 
ment - within the numerical precision - with Vc/N where 
the cusp ofx{v/N) shows up. Through the computation 
of the w-dependence of a topologic invariant, the hypoth- 
esis of a deep connection between topology changes of 
the {El,} and phase transitions has been given a direct 
confirmation. Moreover, we found that a sudden "second 
order variation" of the topology of these hypersurfaces 
is the "suitable" topology change - mentioned at the be- 
ginning of the present Letter - that underlies the phase 
transition of second kind in the lattice model. There 
is no reason why the results presented here should be pe- 
culiar only to the chosen model, and therefore they point 
to a general validity of the relationship between topol- 
ogy and phase transitions, opening a wide field of future 
investigations and applications. 
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FIG. 1. Numerical computation of the Euler characteristic 
for 24 dimensional spheres, v is the squared radius. 

FIG. 2. Id (/p* model. Relative variations of the Euler char- 
acteristic of E„ vs v/N (potential energy density). Lattice of 
iV = 1 X 49 sites. Full line is a guide to the eye. 

FIG. 3. 2d model. Relative variations of the Euler char- 
acteristic of E„ vs v/N (potential energy density). Lattice of 
N = 7 X 7 sites. The vertical dotted line corresponds to the 
phase transition point. Full line is a guide to the eye. 

FIG. 4. 2d (p* model. Relative variations of the Euler char- 
acteristic of E„ vs v/N (potential energy density). Lattice of 
N = 9 X 9 sites. The vertical dotted line corresponds to the 
phase transition point. Full line is a guide to the eye. 
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